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Abstract 

The addition of new multiplets of fermions charged under the Standard Model gauge group is 
investigated, with the aim of identifying a possible dark matter candidate. These fermions are 
charged under SU{2) x U{1), and their quantum numbers are determined by requiring all new 
particles to obtain masses via Yukawa couplings and all triangle anomalies to cancel as in the 
Standard Model; more than one multiplet is required and we refer to such a set of these multiplets 
as a polyplet. For sufficiently large multiplets, the stability of the dark matter candidate is ensured 
by an accidental symmetry; for clarity, however, we introduce a model with a particularly simple 
polyplet structure and stabilize the dark matter by imposing a new discrete symmetry. We then 
' explore the features of this model; constraints from colliders, electroweak precision measurements, 

the dark matter relic density, and direct detection experiments are considered. We find that this 



model can accommodate a viable dark matter candidate. 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics has demonstrated astounding success in 
describing phenomena at energy scales below a few hundred GeV. However, despite this 
success, there are significant questions which the SM does not address adequately. Some 
of these, such as the hierarchy problem and the origin of flavor, arise as theoretical issues 
within the SM, while others, such as the tiny but nonzero value of the cosmological constant, 
the origin of neutrino masses, and the identity of dark matter (DM), reflect discrepancies 
between observation and the SM prediction. All of these questions indicate that, despite 
the remarkable success of the SM, it is not a complete theory of nature. 

Given both the successes and the shortcomings of the SM, it makes sense to consider 
models of new physics which address one or more of the above questions while not radically 
altering the structure of the SM. Following this philosophy, we investigate the possibility 
of adding new fermions, including a possible DM candidate, to the SM without altering its 
gauge group, SU{3) x SU{2) x U(l). (For examples of works where exotic fermion multiplets 
have been added to the SM, see While all fermions in the SM transform as singlets or 

as the fundamental representations under the SU (3) color and the SU{2) weak interactions, 
higher-dimensional representations are possible. Here, we investigate the possibility of new 
fermions which transform as higher- dimensional multiplets under SU{2). We will, in fact, 
add more than one new multiplet to the SM, and we will refer to this set of multiplets as a 
"polyplet" . 

There are several advantages to adding higher- dimensional SU (2) fermion multiplets to 
the SM. The first is that we can hope to arrange the quantum numbers of the multiplets 
such that they obtain their masses only via Yukawa couplings with the SM Higgs boson; 
thus our model would contain no new mass scales and would predict all of its new particles 
to be near the electroweak scale. Second, if the dimensionality of these new multiplets is 
large, with isospin greater than 1, it is impossible to write down dimension- four mass terms 
which mix these new fermions with those of the SM. In addition to simplifying the model, 
this implies that the lightest member of the polyplet cannot decay to SM particles through 
renormalizable operators; this is useful as we hope that our polyplet will contain a DM 
candidate^ For related works in which SU{2) multiplets have been added to the SM in 
order to solve the DM puzzle, see jsj-joj. 

Several considerations must be taken into account when introducing new SU{2) multiplets 
to the SM. First, constraints on the existence of new neutral or charged massless particles 



from colliders [10[ and cosmology [11|, |12| imply that all members of these new multiplets 
must be given masses. We choose the quantum numbers of these new multiplets such that all 
masses are acquired via Yukawa couplings with the SM Higgs; we wish to avoid dimension-3 
mass terms and we disallow mixing with the SM fermions. Additionally, we wish to avoid 
the possibility of a charged stable particle (see and references therein for constraints 



on charged relics), so we insist that the lightest new particle is electrically neutral. Finally, 
as we are interested in possible DM candidates, we only consider colorless particles. 

After giving masses to our new fermions, we explore their phenomenology. Because 
their transformations under the SM gauge group are specified, their interactions are largely 
predetermined, although some couplings are tunable via mass mixing. We consider present 
experimental bounds on these new higher- dimensional multiplets and find that they are 

^ We note, however, that, for simphcity, we will work with a polyplet whose largest isospin is 1 and introduce 
a new conserved quantum number to disallow the unwanted operators. 
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allowed for a wide range of new particle masses. We also find that our DM candidate can 
achieve the observed relic density via annihilation to gauge bosons and top quarks while 
simultaneously evading direct detection constraints. 

We organize the remainder of this paper as follows. In Section [ITl we introduce our model, 
and discuss its particle content, couplings, and masses, and restrictions put on the model by 
requiring anomaly cancellation. We then consider constraints from searches at colliders and 
precision electroweak measurements in Section lllli In Section IIVI we investigate constraints 
placed by DM direct search experiments, and, finally, in Section |Vl we conclude. 



II. THE MODEL 

A. Notation and Basic Considerations 

We begin this section by specifying our notation and discussing some basic considerations 
which must be taken into account in order to make all of the particles in our model massive. 
The polyplet we will consider contains four fermion multiplets, charged under SU{2) x f/(l), 
which we will denote as left-handed chiral spinors Xni,Yi- Here we label our multiplets with 
their number of components (where rii = 2/j + 1, with Jj the multiplet isospin) and 
hyper charge Yi. 

In order to evade the constraints on the existence of new massless particles, we must 
arrange masses for all of the particles in our model. At the same time, we wish for our 
model to be predictive; therefore, we insist that all fermion mass terms be dimension-4 
Yukawa couplings to the (hypercharge 1/2) SM Higgs doublet H 

XnuY, X Xn„-y,-i X Hi, 

Xn,,Y, X Xn,,-y,+i X H\, (l) 

where x is the multiplication needed to render the mass terms SU (2)-invariant and the value 
of Yj has been determined in terms of Yi by hypercharge conservation. After electroweak 
symmetry breaking, 

where v ~ 246 GeV; we thus expect the fermion masses to be tied to the electroweak scale. 

Because the SM Higgs field is an isospin doublet, two x multiplets can couple to each 
other via a Yukawa term only if their isospins differ by ±1/2 or, equivalently, if nj = nj ± 1 
in Eq. ([T]). Thus, x multiplets with half-integer isospin (or even rii) will only be coupled 
to those with integer isospin (or odd n^), and vice- versa. This implies that, in order for all 
particles in these multiplets to obtain masses, we must have an equal number of particles in 
even-nj multiplets and in odd-n^ multiplets. 

Let us consider the possible Yukawa couplings of a multiplet Xn,Y with n components. 
Let us assume it couples to another multiplet Xn-ix'- 

Xn,Y X Xn-l,Y' X H^*l. (3) 
^ 2 

However, as Xn,Y contains one more component than Xn-iy-, this cannot render every 
component of Xn,Y massive. Thus, Xn,Y must couple to another multiplet. We can either 
couple Xn Y to another multiplet with n — 1 components, or to one with n + 1 components. If 
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we couple it to another multiplet with n — 1 components, our fourth multiplet must contain 
n — 2 components in order to balance the number of particles in even-n and odd-n multiplets, 
and our model contains multiplets of size n, n — 1, n — 1 and n — 2. If we instead couple Xn,Y 
to a multiplet with n + 1 components, we obtain multiplets of size n, n, n — 1, and n + 1. 
These two choices are clearly degenerate by a redefinition of n. We therefore can write the 
four members of our polyplet as: 

(4) 

Next, we note that, for a given value of Y~^, the members of Xn+i,Y+ will have n + 1 
distinct values for their electric charge, Q = /(s) + Y. In order for all of these fields to 
become massive via Yukawa couplings with Xn,Y and Xn.y, these same values of Q, but with 
opposite sign, must be represented in the components of Xn,Y and Xn,Y' ' which requires that 
Y 7^ Y'. We therefore take the Yukawa couplings of Xn+i,Y+ 

Xn,Y X Xn+1.Y+ X -f^* i, 
2 

Xn,Y' X Xn+1,Y+ X Hi, (5) 

from which we can deduce that 

r + r+ - i = 0, 

F' + F+ + i = 0, (6) 
Y' = Y -1. 

Now we consider the Yukawa couplings involving Xn-i,Y~- There are three possibilities 
for these couplings. First, Xn-i,Y- could couple to both Xn,Y and Xn,Y'- ^^i^ case, we 
yield couplings analogous to those in Eq. ([5]), and we trivially obtain Y~ = Y~^^ However, 
it is also possible that Xn-i,Y- couples only to Xn,Y or only to Xn,Y'- In these cases, one can 
have either a coupling of the form 

Xn,Y X Xn-1,Y- >< Hi, (7) 

or one of the form 

Xn,Y' X Xn-1,Y- ^ H*_i, (8) 

2 

which yields Y~ = Y~^ — 1 or Y^ = Y^ + 1, respectively. We will thus adopt the following 
notation for our fermion multiplets: 

Xn,y, Xn,y-1, X„+l,i-y, Xn-1,^-Y+A^ (9) 

where A = 0, 1 or —1. 

With these preliminaries out of the way, we will now explore the constraints on n, and 
the multiplet hypercharges imposed by anomaly cancellation. 



^ It is also possible to include only one of these two couplings, but this will not change the resulting 
relationships between the multiplet hypercharges. 
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B. Anomaly Cancellation 



When building models which contain new particles or interactions, one must take care 
not to ruin the cancellation of triangle anomalies which occurs in the SM. (For a related 
work on anomaly-free sets of fermions, see [l6|.) As our multiplets are not charged under 
SU{3), there are only three anomaly relations relevant here: f/(l)-gravity-gravity, U{1)^, 
and U{1) X SU{2)'^. These three relations take the form 

J2 ^^^^ = 0' 

multiplets i 
multiplets i 

E ^(^)^^ = 0' 

multiplets i 

where C{i) = \k{k + l)(2/» + 1). 

When we apply these constraints to the quantum numbers of our multiplets shown in Eq. 
OH]), we obtain 

(n-l)A = 0, 

1 n 1 

3(n - l)A{- - Yf + 3(-- + A\n - 1))(- -Y) + {n- l)A^ = 0, (11) 

(/(n + 1) + f{n - 1) - 2/(n))(^ -Y) + f{n - l)A = 0, 

where, for clarity, we have written the third relation in terms of f{x) = (x — l)x{x + 1)/12. 
Note that /(n + 1) + /(n - 1) - 2f{n) = n/2. 

Let us examine the first relation in Eqs. (fTTj) . This equation has two solutions; n = 1 
or A = 0. Plugging either of these solutions into the other two equations in (fTTl) . we 
obtain Y = 1/2. However, we note that, for Y = 1/2, odd n implies that our multiplets 
contain particles of half-integer charge. As such fractionally-charged (and stable) particles 



are observationally ruled out |l3Hl5|. we insist on even n and, therefore, A = 0, yielding the 
multiplet^ 

' Xn+1,05 Xn-1,0- (12) 

This set of quantum numbers, however, allows dimension-3 mass terms like mxn-ifl x 
Xri-1,0) ^Xn+ifi X Xn+1,0) and mXn +i X Xn -i) where again x is the multiplication needed to 
make an SU{2) singlet. We would like to disallow these terms, in order to have a predictive 
model with masses set by the electroweak scale. We can do this either by introducing 
additional fields to loosen the constraints set by anomaly cancellation, or by assigning a new 
conserved charge to our multiplets. For simplicity, here, we will impose the conservation 
of a new charge, which we refer to as "pseudo-lepton number". A discussion of how we 
can instead loosen the restrictions imposed by anomaly cancellation via the introduction of 
additional fields is given in the Appendix. 



^ We note that for even n, Tr lf„s, when summed over aU multiplets, yields an integer; thus, this choice also 



satisfies the Witten anomaly condition |17| . 
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With all of the hypercharges now specified, we briefly discuss whether or not the four 
allowed Yukawa couplings are sufficient to render every member of the polyplet massive. To 
see this, we revisit the couplings of Xn+i,Y+ given in Eq. ([5]), which now take the form 

y X„,-l X Xn+lfi X Hi, (13) 

and 

y' Xn i X Xn+1,0 xH*_^, (14) 

where we have introduced two Yukawa couplings y and y'. The couplings in Eq. ( !T3|) will 
couple the components in the first n components of Xn+i,o, 

CiXn,_i{n)xn+i,o{l) + C2Xn,-iin " l)Xn+i,o (2) . . . . + C„x„,_ i (l)x„+i,o (^) , (15) 
while those in Eq. (IT^ will couple the components of x„,i the last n components of Xn+i,o, 

C2Xn,i('^)Xn+l,o(2) + C'^Xnli^ " l)Xn+l,0 (3) • • . . + C^^iXn,! (l)Xn+l,o(^ + 1), (16) 

where the coefficients Cj and C- are products of f/v^, Clebsch-Gordon coefficients, and 
either y or y', respectively. Two members of Xn+i,o, Xn+i,o(l) and Xn+i,o{n + 1), couple to 
only one component of either Xn ^ Xn-^ ^"^^ thus mass eigenstates. However, for 
1 < i < n + 1, o(i) couples to a linear combination of one component from each n-plet, 

(Ca„,-i(^ + 1 - + C:x„ i(^ + 2 - z)) X„+i,o(^). (17) 

while the orthogonal combination of Xn -i(^ + 1 ~ and Xn i('^ + 2 — i) remains massless. 
When we include the analogous couplings of Xn-ifl, 

(AXn,-i in + l-t) + D[xn,l{n + 2 - «)) Xn-i,o(^ - 1), (18) 

these previously massless combinations will receive Yukawa couplings as long as D[/Di ^ 
C'JCi. As the coefficients Cj, Di, D[ are functions of arbitrary yukawa couplings, it is 
reasonable to expect that they can generally be adjusted to render all of the fields massive, 
as long as none of the relevant Clebsch-Gordon coefficients are zero. 

Up until this point, we have not specified n. We will take the simplest polyplet which 
yields particles with integer electric charge, n = 2. We assign X2,+J- X2,-i ^ pseudo- 
lepton number of 1 and X3,o and xi,o a pseudo-lepton number of —1. All SM particles have 
a pseudo-lepton number of 0. We note that these assignments also prevent our multiplets 
from mixing with the SM leptons. 

As a side note, we briefly mention that the quantum numbers for the SM quarks satisfy 
the relations given in Eq. for n = 1, Y = 1/3; here, Xi i; Xi and X2 1 correspond to 

I — I '3 '3 '6 

d^, M^, and the SU{2) quark doublet, respectivelyo Additionally, the SM leptons form an 
incomplete n = 1, Y = 1 polyplet, with xi,i corresponding to e^, and X2,-i to the lepton 
doublet; the addition of xi,o (i-e., a right-handed neutrino), would give the neutrino a Dirac 
mass. Of course, in the SM, anomaly cancellation is achieved by having contributions from 
the quark and lepton sectors cancel against each other, unlike in our case here. 

We now move on to discuss some of the phenomenology of the n = 2 model introduced 
above. 



Note that for n — 1, there is no fourth muhiplet, as it would have n — 1 = components. 
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C. Features of the model 



After the dust has settled, our model contains 8 chiral states. X2,+i X2,-i ^ach 
contain one charged and one neutral state, 



2 '2 




= , (19) 



where we have introduced superscripts to indicate the electric charge of each member of the 
multiplet. Being SU{2) doublets, X2,+i -^2,-5 ^^^^ couple to the W and Z bosons. On 
the other hand, xi,o (= Xio); being uncharged under both SU{2) and f/(l), will not couple 
to the 7, Z, or W. Lastly, X3,o contains two charged and one neutral state, 

X3,o = Xlo ■ (20) 

V ^3,0 / 

While all three states of X3,o will couple to the W, we note that Xs.o; having both J(3) and 
hypercharge 0, will couple to neither the 7 nor the Z. 

Next, we explore the Yukawa couplings of these fermions. First, we write down the 
couplings which involve xi.o, 

^ ^ -y+i (x2,+i e H xi,o) - y~i (x2,-i e H xi,o) • (21) 

where we have labelled the Yukawa couplings yij according to the sign of the hypercharge of 
the doublet X2±i subscript 1 refers to the number of components in the hypercharge- 

multiplet. H denotes eH* where e = iT2 and T2 is a Pauli matrix. After electroweak 
symmetry breaking, — )■ ^ ^ ^ and — )■ ^ ^ ^ ^ , and we obtain the mass terms 

-(-l)^^2,+iXi,o - ^X2,_iXi.o- (22) 

The Yukawa couplings involving the isospin triplet Xa.o are slightly more complicated. We 
need terms of the form 

- y+3 X2,+l X X X3,0, 

-y-3 X2,-i X H X X3,o, (23) 
which form isospin singlets. These give mass term^ 

^ ^ (^J+I^^o ~ ^4,+iX3,o) - ^ (-^4-1X3,0 + Xl_iXto^ ■ (24) 
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The form of these mass terms is simple to see using the following considerations. Our mass terms contain 



two isodoublets and one isospin triplet. Two isospin doublets X — ( '^^ | ^^J^d x' = ( ''''^ | can be 

combined into the triplet object xtXi+XiX^ _ rp-j^^ Clebsch- Gordon coefficients to combine two 

\ XiXl / 

SU{2) triplets into a singlet are -i^ , — , -i^ . Using these relations, normalizing the Yukawa couplings 
to remove the factor of and setting H and H i^o their vacuum expectation values gives Eq. (|24|) . 



Thus, our model contains four massive particles. We have two charged particles, with masses 
^^^^ and ^^^1^, respectively, and two neutral states, whose masses are determined by the 
mass matrix M in 




A priori, the Yukawa couplings in Eq. fl2^ are arbitrary complex numbers. However, we 
note that we can remove three of the phases in these Yukawa couplings via U{1) rotations 
on the multiplets. To see this, we write Eq. (125|) in terms of the mass eigenstates, which 
are related to the neutral fields in the polyplet via unitary transformations. We write 



{mN,) = i^ Xl. <_| j V\ (26) 

and 

\ _rT( Xi.o 



where U and V are unitary matrices. We will write U in the form 

where |ap + |/3p = l. Via a redefinition of 6u, we will in turn write this as 

^-""(L-l.)(-^")(oA)- 

where su and cu are real and + cIj = 1. Similarly, we will express V as 

v'--(jA.)(-x:)(JAO- (^o' 

Then, writing the mass terms as 

M in Eq. fl25|) becomes 

^^../l \ /cv^ -sv^\ /I \ /mi \ / /I \ 

\^ e^-^^^ J 1^ cv J \ e'" J \ m2 J \ -su cu J \0 e~''^v^ ) ' ^ ^ 

where we have absorbed 9u and into a single phase 9 and similarly combined the Oyi and 
6'[/i phases into a single angle a. By performing field redefinitions on our multiplets, we can 
absorb the e*^, e^'^^^, and e^'^^'^ factors into the fields. (Note that we do the field redefinition 
on the entire multiplet, and not just on the neutral component, to avoid introducing unnec- 
essary phases in charged-current couplings.) Thus, the mass matrix can be parameterized 
by five quantities, mi, m2, cy, C[/, and a. 
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For our purposes, however, it is more useful to retain e*^, e 



^<PV2 



and e^'^^'^ and use them to 



write the mass matrix in a more convenient form. In order to have standard mass terms for 
the charged states in Eq. ( 12^ . we wish to have 1/+3 and ?/_3 real and positive. Additionally, 
for convenience, we will define our fields such that y+i is real and negative, and hence the only 
complex phase in the Yukawa matrix in Eq. (1251) is contained in y_i. This is accomplished 
by keeping M of the form in Eq. ( 15^ . but now 9, (j)v2 and (f)u2 are functions of the other 
five parameters; they will fall out of couplings and thus will not contribute to any of the 
physical quantities discussed below. 

Now that we have a general form for our mass matrix, we wish to verify that we can 
choose Yukawa couplings such that the lightest particle in our polyplet is neutral in order 
to have a possibly viable DM candidate. To see this, let us take the simple case —y+i = 



y+3 



/V2 = y^s/V2 = 1, 



y-i 



e**^, where is a phase. In this case, the two charged states 



each have mass v, while the neutral states have mass squared of ± 
For (f) = 71, this gives degenerate neutral particles with mass v. However, for 



1 — 2(1 — cos 



0, we 



obtain mi = 0, m2 = \/2v. By perturbing about this latter solution, we can make the 
lighter of the neutral states massive yet still lighter than both charged states. Of course, one 
can also consider more general Yukawa couplings than those given in this simple example. 

We now consider the couplings of our particles to the SM gauge bosons. We first note 
that the Z does not couple to xi,o or X3,o both have Y = 7(3) = 0. The Z does, however, 
couple to the charged states, as well as the neutral components of X2 i and X2 - i- 



COSOw '2 "^'2 \2 J ^'2 ^'2 \^ 



^'2 ^'2 
t; 



-- + sm 



w 



+xtl(r''xto (1 - sin^ Ow) + X3>^X3,o (-1 + sin' Ow)] 



(33) 



where 9w is the weak mixing angle, a'^ = (1, — and "0^ is the set of Pauli matrices. We 
note that the couplings of the neutral fields are equal in magnitude but opposite in sign; we 
will see below that this allows the coupling of the Z to the neutral mass eigenstates to be 
somewhat tunable and allows the Z to couple to states which are off-diagonal in the mass 
basis. 

The charged- current interactions of the doublets X2 ^ X2-^ tak.e a form similar to 
that of the left-handed fermion doublets in the SM, 



£d4 
V2 



W;^ixt\^'xU+xl^ ^) + W-ixl\crX^+X,^ rcr'^xl.,) 



(34) 



while our SU (2) triplet field interacts via 



W. 



, {xto^'X3,o + Xtyxlo) + W'^ 



iXsl^'xlo + Xto^'xto) 



(35) 



We now re-express Eqs. (13311351) in terms of the mass eigenstates. For the neutral states, 
we combine A^i and A'^^ into a single massive Dirac spinor with Ni as its left-handed compo- 
nent. This Dirac spinor has mass mi and we will now refer to it as Ni. Analogous notation 
applies to A^2- Without loss of generality, we will assume that 777.2 > ^i- The couplings of 
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these particles to the Z boson are 



cos 6' 



w 



(36) 



where Pi^ is the left-handed projection operator. Here we see that the mass-diagonal cou- 
plings of Nx and N2 to the Z can be tuned to be small if \cy\ ~ |sy| and that the off-mass- 
diagonal coupling can be substantial. 

Similarly, the two massive charged particles also interact with the Z . Denoting the 
particle with mass as Xi and the particle with mass as X2, their couplings to the 
Z are 



cos Q ^ 



+X2r 



Q - sin^ Ow^ Pl + (1 - sin^ Ow) Pr 
-i + sin2 ewj Pl + (-1 + sin^ Ow) Pr 



X2\ 



(37) 



Note that, with our definitions, Xi is positively charged and X2 is negatively charged. 
In the mass basis, the charged-current couplings in Eqs. and become 

C D -^W; [xi1^Pl{cvNx - sye-^"iV2) + {syNi + cve'''N2)l^PLX2] + h.c, (38) 



and 



C D gW; [-Xi1^Pr{suNx + c^/iVa) - {suN^ + Cf/iV2)7^Pi?X2] + h.c. 



(39) 



respectively. 

Now that we have established the couplings of the particles in our polyplet, we can discuss 
some of the constraints on this model. 



III. CONSTRAINTS 



As this model contains new weak-scale fermions which couple to the and Z bosons, 
it is important to examine existing constraints from colliders. A basic constraint comes 
from the measurement of the invisible width of the Z at LEP [isj, which requires that a 
new stable neutrino with SM couplings to the Z have a mass of at least 45.0 GeV. We 
thus expect that our lightest neutral particle would have a mass of at least this magnitude, 
although this constraint could be evaded if only one neutral particle was light and |sy — Cy| 
was sufficiently small. 

There also exist a few constraints from colliders on new neutral particles with masses 
greater than The authors of 19j use the single-photon spectrum from DELPHI 

2ol . 21 1 to constrain the production of invisible particles in e^e~ collisions via either effective 



operators or on-shell mediators. For the case most relevant here, production of Xi of mass 
> Mzjl via a Z boson, however, the limits are weaker than those from direct detection 
which will be discussed in Sec. HVl Similar conclusions apply to studies of mono-jet events 
at the Tevatron [22,] . 
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It is also possible that the heavier of our two neutral states, A^2, could be produced at 
colliders via Z — )■ N2N1 or Z — )■ N2N2- If the off-diagonal Z couplings are not tiny, N2 
could decay down to A''i via N2 — )■ NiZ. (Also, if N2 is heavier than xi oi' X2, it could 
decay to Xi,2W^^, followed by the subsequent decay of the Xi,2-) The similar decay of fourth- 
generation neutrinos with both Dirac and Majorana masses has been investigated in j23| . 
which used LEP SUSY searches jl^ to conclude that LEP would likely exclude masses of 
the heavier neutrino below ~ 100 GeV, and possibly as high as 130 GeV, if the lighter of 
the two states had a mass of ~ 50 GeV. Although the couplings in our case are somewhat 
different, we would expect the limits here to be roughly similar. For a related analysis at 
CDF which gives weak constraints, see (25| . 

There also exist collider limits on heavy charged leptons. L3 [26| considered the case of 
a new charged lepton which decays to a heavy neutral lepton L° and a W^, assuming 
that the L° weighs at least 40 GeV and that the mass difference between the and the 
LP was between 5 and 60 GeV; they place lower bounds on the mass of ~ 100 GeV. As 
in the case of heavy neutral leptons, the constraints that these limits place on our model 
are dependent on the mixing angles, but we would nonetheless expect these constraints to 
roughly hold. 



One may worry that this model would be ruled out by recent Higgs searches at LHC [27 



28| or at the Tevatron [29[ in the context of fourth-generation models, as models with extra 
fermion generations often have signals which are boosted relative to the SM, particularly 
in the case of if — ?• 77, which is mediated in part by fermion loops. (For a recent survey 
of the experimental limits on such scenarios, see |30|.) However, this model contains no 
heavy colored particles beyond the SM top quark; therefore, there is no enhancement in 
Higgs production via gluon fusion in this model. Additionally, while this model does give 
a contribution to the decay H — )■ 77, this contribution, being from a fermion loop, is 
opposite in sign to the main SM contribution, which is mediated by a loop. Therefore, 
the branching fraction for — t- 77 is actually reduced in this model relative to the SM. 
Additionally, if the Higgs can decay to non-SM particles, then this can loosen the current 
Hi ggs exclusion bounds, similar to the case of Higgs decays to fourth-generation neutrinos 
in |3l|. 

In addition to constraints from direct searches, however, we may consider the implica- 
tions of heavy fermions on electroweak precision parameters. Thus, here we consider the 
contributions of our new fermions to the oblique parameters 5* and T [s?]. While a full 
exploration of the parameter space is beyond the scope of this work, we consider a special 
case, with a mass matrix of the form 




1 / —mi —1712 

y/2\ mi -1712 



(40) 



which gives A^i a mass of mi and N2, Xi; ^ind X2 a mass of m2. Diagonalizing this mass 
matrix leads to the mixing parameters 



(41) 
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FIG. 1. The contribution to S in the simpHfied mass and mixing scenario as a function of nii for 
three values of the mass difference m2 — mi. 

Inspecting Eqs.f l36p . we see that this scenario yields no mass-diagonal Z couplings to Ni 
and A^2- This feature will be particularly useful when we consider the potential of A^^i to be 
a dark matter candidate in Sec. HVl Additionally, from Eqs. fl38|) and fl39|) . we see that the 
charged-current couplings of A^i are entirely left-handed. 

In Fig.([T]), we plot S* as a function of mi for this simplified scenario of masses and mixings. 
Here, we take = mi + 10, 100, 200 GeV. We see that S is typically C(0.01 - 0.1) in all 
three cases. We note that in the case that mi = m2 >> Mz, S = 1/ (Sti). We compare this 
with the results of 33 ^ which uses results from the Gfitter Group [Isj to find that the new 
physics contributions to S and T are constrained to be 



^1 



T\ 



u=o 



0.07 ±0.09 



In- 



mt 



u=o 



0.10 ±0.08 + 



127r (120GeV)2 
3 



In 



mt 



167r cos2 9w (120GeV)2 ■ ^^^^ 

Additionally, for general mi, m2 in this simplified scenario, T = 0; this is because the Yukawa 
couplings do not break the custodial symmetry for the case y+i = y_i, ?/+3 = y-sIZl (Outside 
of this specific scenario, however, T is not generally zero.) We thus see that this model can 



^ Also see [s^l for a discussion of the oblique parameters in the context of fourth-generation scenarios. 



^ We can see this by rewriting our fields as X3,o — 




X2 



H 



^ 2 2 ' 

Under SU{2)l and SU{2)r rotations Ul and Ur, X3,a ^ UlX?.,qUI, X2 ^ UlX^U]^, 



H — > UlHU\i. Yukawa couplings in this scenario can be written as Tr[x3^oX2-ff^] and Tr[x2^^^]xi.o; when 
the H vacuum expectation value breaks SU{2)l x SU{2)ii down to the custodial symmetry, Ul = Ur, 
these terms remain invariant. 12 



easily produce values for 5* and T which are in reasonable agreement with experiment for 
low Higgs masses. We will return to this in Sec. IIVI when we discuss the interplay between 
these results and constraints from direct detection. 

Lastly, we briefly mention a few other possible constraints on the polyplet model. First, 
as the couplings of our new fields contain a complex phase a, it is possible that our model 
contains contributions to CP- violating observables, such as electric dipole moments (EDMs) 
of SM particles (constraints on which can be found in [iq|), once we move beyond tree- level 
calculations. We note, however, that a can be chosen to be small; we thus do not consider 
constraints from EDMs here, but they could be of use in a more specific model. Also, 
as has been pointed out in |36|, the addition of new heavy fermions to the SM can have 
ramifications for vacuum stability; here, we simply assume that this issue will be resolved 
by some other new physics added to the SM. 

We now move on to discuss the viability of Ni as a DM candidate. 



IV. Ni AS DARK MATTER 

As the lightest member of our polyplet, A^i is stable and therefore a potential dark 
matter candidate. (For large n, this would be a result of there being no dimension-four 
S'?7(2)-invariant operators which couple the members of our polyplet to the SM fermions; 
for n = 2, it is enforced by our choice of pseudo-lepton number.) In order for A^^i to possibly 
be a realistic DM candidate, it must fulfill two basic requirements. First, it must annihilate 
to SM particles in the early universe sufficiently to have a relic density which is equal to 
or less than the observed DM density today, and, second, it must interact with nucleons 
sufficiently weakly that it evades current constraints from direct detection. Both of these 
processes are mediated via interactions of A^i with the Z boson, the Higgs boson, or a 
combination of the two. 

In trying to simultaneously fulfill the relic density and direct detection requirements, a 
potential conflict arises. If DM is a thermal relic, its observed density of VtoMh? = 0.112 ± 
0.006 implies a thermally-averaged annihilation cross-section of < (j\v\ >~ 3 x 10~^® cm^/s. 
(For a general review of the status of dark matter and its detection see, for example, 10|.) 
The observation that this cross-section is typical of weak-scale processes is commonly referred 
to as the WIMP miracle. However, direct detection experiments have put constraints on 



DM-nucleon cross- sect ions as stringent as C(10-^^cm2) ~ 0{2 x 10"^^ GeV'^) [37|, [38 
If DM-nucleon interactions are mediated b y p hysics at a new scale A, DM-nucleon cross- 
sections would typically be of order /i^/A"^ |39|, where is the DM-nucleon reduced mass, 
0{1 GeV). From this, one can see that a value of A ~ 100 GeV or even A ~ 1 TeV would 
produce a DM-nucleon cross-section which is ruled out by direct-detection experiments by 
several orders of magnitude. Thus, in order to have a DM candidate which simultaneously 
has an acceptable relic density and evades the direct detection bounds, there must be some 
signiflcant difference between the physics which controls the DM annihilation and that which 
dominates DM scattering with nucleons. 

As a flrst example, we consider the case of a relatively light Ni, mi < mv^, where the A^i 
can only annihilate to the light SM fermions. As long not very small (and 

annihilation via a if is non-resonant), annihilation via a Z will dominate over if- mediated 
processes. If we take the Z coupling to be full-strength, |sy — Cy| ~ 1, Ni will comprise a 
subdominant component of DM, yet will fail to evade direct detection constraints. Also, as 
long as both DM annihilation and nucleon scattering are dominated by diagrams mediated 
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by a Z boson, adjusting the coupling of A^i to the Z will not help mitigate this conflict. 
Although reducing Sy — Cy will reduce the effective coupling between the DM and nucleons, 
it will similarly reduce the annihilation cross-section, thus increasing the current-day relic 
density; this effect will compensate for the reduced coupling, leaving the direct-detection 
constraints largely unchanged. Therefore, we must have non-negligible contributions to A^i 
annihilation from processes mediated by the H. 

If we keep mi < M-\y but decrease the mass-diagonal coupling between the Z and A^^i to 
zero such as in the simplifled scenario we explored in the discussion on oblique parameters 
in Sec. IIIIl we can have A^^i interact with the SM entirely via the Higgs. In this case, one 
might hope that the direct-detection cross-section will be sufficiently suppressed by light 
quark Yukawa couplings to evade the relevant constraints. If mi is small, however, A^^i can 
only annihilate to the light SM fermions, and, thus, the dominant piece of this cross-section 
is suppressed by the Yukawa coupling of the b quark. Unless mi ~ mu/'^i this annihilation 
cross-section is too small to yield a reasonable relic density. Additionally, if one chooses mi 
to obtain the correct relic density, direct detection bounds still rule out this scenario for a 
light A^i. The direct detection cross-section is dominated by the s-quark contribution to the 
nucleon, 

2 / \ 2 

(^nudeon-Ni ~ -^TVs 1 ~o 9~ 5 ('*"^J 

where Bfq^ = 3.36 ± 1.45 is the factor which converts between quark-level and nucleon-level 
effective operators; we take this value from 22|] which used m^ = 105 ± 25 MeV. (We note 
that the error bar on this factor is large and different conversion factors are used in the 
literature; see, for example, |40i, 41]. Additionally, the expected contributions from m, d, and 
heavy quarks are smaller than or comparable to the error on the value in Eq. (I43l) .) This 
yields 

« 3.7 X 10-« cm= (^^) (— — j (44) 

with an error bar of the same order. As the limits on the spin-independent DM cross-section 
from XenonlOO 37j are (9(10~^^ cm^), this scenario is ruled out by approximately an order 
of magnitude for low Higgs masses. 

If we are willing to allow A'^i to comprise only a very subdominant component of DM, we 
can use resonant annihilation via the Higgs to make the current relic density of A'^i small. In 
this case, we tune the A''i mass to be mi ~ Mh/2, but we do not flne-tune the Z coupling to 
be small. Although the A'^i-nucleon cross-section will not be suppressed, one can hope that 
the reduced Ni relic density will allow A'^i to evade the direct detection constraints at the 
expense of identifying the primary component (s) of DM. In Fig. [21 we show the relic density 
Qisfih"^ of Ni as a function of mi assuming resonant annihilation through H for mn = 115 
GeV; we include the effects of thermal averaging similarly to what was done in (See 



also |43l. |44| . ) (We note, however, that only for points very close to the minimum does the 
resonant contribution to the cross-section dominate over annihilation through a Z, which 
is not included in the plot.) We see that it is possible to achieve relic densities as low as 
few X 10""^. However, this suppression is inadequate to evade direct detection constraints; in 
Fig. Owe plot the A''i-nucleon cross-section calculated for XenonlOO weighted by Qiy-^^/QDM 
and flnd that it is still ~ 2 orders of magnitude too large to be compatible with the results 
of XenonlOO [sj. 

A solution is possible, however, for larger mi. If we tune the mass-diagonal coupling of A^i 
to the Z to zero {sy — Cy must be < 5 x 10"^ to evade direct detection limits), and increase 
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FIG. 2. Qj\i-^h'^ as a function of mi in the case of resonant annihilation through H. rriH = 115 
GeV. Effects of annihilation via a Z are not included and dominate away from the minimum. 
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FIG. 3. The A^i-nucleon cross-section weighted by Q7Vi/^^dj\/) calculated for XenonlOO, in the case 
of resonant annihilation to light fermions via H. mn = 115 GeV. 

mi so that annihilations to gauge bosons NiNi — )■ if — )■ W'^W" , ZZ and ti become possible, 
we can simultaneously fulfill the relic density requirement and evade the constraints from 
direct detection. In this case, the annihilation cross-section for Ni is not suppressed by any 
small Yukawa couplings, while the main contribution to the direct-detection cross-section is 
still suppressed by the s quark Yukawa. This solution does, however, require a rather heavy 
Higgs boson, niH ~ 0(600 GeV). We plot the relic density of Ni in this case in Fig. HI for 
niH = 600, 700 and 800 GeV, along with the observed DM relic density VtuMh'^ ~ 0.11. 
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FIG. 4. A^i relic density as a function of mi for rriH = 600,700 and 800 GeV. The horizontal line 
is ^DWih'^ ~ 0.11. 



For these values of nii and itlh-, it can be seen from Eq. (jHj) that typical values of 
the direct detection cross-section are (9(10~^^ cm^). Thus, we see that if mi ~ 200 GeV, 

— Cy < 5 X 10~^ and we have a heavy Higgs, A^^i can be a successful thermal DM 
candidate; additionally, the typical A^i-nucleon cross-sections are accessible by near-future 
direct detection experiments. 

We briefly comment on the feasibility of a heavy Higgs boson. Values of itlh up to 
600 GeV have been ruled out by CMS [23|, although it is possible that this constraint 
could be loosened slightly in the polyplet scenario as new decay modes of the Higgs, H — )■ 
NiNi, N2N2, X2X2 are available. A more concerning constraint comes from the oblique 
parameters. If we again consider Eq. (l42l) . we see that a 600-800 GeV Higgs boson requires 
that new physics make a contribution to T of approximately 0.35-0.4. To see that such a 
contribution can arise in the polyplet scenario, we revisit our simplified couplings in Eq. 
f l4ip . but we now take Su to be small but nonzero. In this case, the masses of the neutral 
particles are still mi and m2, but the charged particles now have masses m2Cu ± rriiSu- A 
plot of T for Su = 0.2 is shown in Fig. [Sj we see that values T ~ 0.4 can be achieved for 
reasonably small su- 

As our solution to the dark matter relic density and evasion of the direct detection 
constraints relied on tuning the NiNiZ coupling to be close to 0, one may worry that 1-loop 
corrections might ruin this solution. Here, we assume that the fine-tuning of the NiNiZ 
coupling includes such corrections. However, other 1-loop effects could also be relevant. It 
is possible that Ni could interact with the SM via a magnetic moment. We have checked 
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FIG. 5. T as a function of mi for su = 0.2 and three values of m2 — mi. 



(using the results of t45i]), for the simpUfied coupUngs in Eq. fj4T|) . that the magnetic moment 
of A^i vanishes at one-loop order, due to a cancellation between diagrams containing a xi 
and those containing a X2- However, we note that the sizes of these individual contributions 

and, therefore, the 



are of the order of those which would be ruled out by XenonlOO [3j 
ability to evade this bound will be somewhat dependent upon the particular mass and mixing 
parameters chosen for the model. It is also possible that A^^i could interact with quarks via 
box diagrams; we assume that these contributions to the direct-detection cross-section are 
suppressed by (9((a/47r)^) relative to a weak-scale cross-section. 

Given that our set of parameters for Ni to be a DM candidate involved a heavy Higgs 
boson, we briefly comment on the scenario in which the current hints from LHC 128| 
and the Tevatron [29| of a Higgs with a mass of ~ 125 GeV are realized. If it is eventually 
confirmed that the Higgs sector consists of a single, comparatively light scalar with SM 
couplings, then constraints from direct detection may require us to consider somewhat more 
complicated solutions. We outline a few of these here. 

If we are willing to accept additional fine-tuning, we can return to the scenario of resonant 
annihilation; from Fig. [3l we can see that if the mass-diagonal coupling of Ni to the Z 
were reduced from its maximal value by approximately one order of magnitude or more, 
constraints from direct detection experiments could be evaded for mi ^ mn/^. 

Alternatively, we could accomplish much the same situation using coannihilations; if the 
mass-diagonal Z coupling were made small, but the Ni was close in mass to one of the other 
polyplet fields (with a mass difference of order a few GeV or less), we could have the relic 
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density controlled by annihilations via the Z and/or the W without the Z contributing sig- 
nificantly to direct detectionjl In this case, A^^i would have to be a subdominant component 
of DM, in order to evade direct detection constraints from its interaction via a 125-GeV 
Higgs (see Eq. (HI])). We note that the constraints in 23, 25, 2^ disappear in the case that 
the mass splitting between A^i and the next lightest state is small. 

Lastly, another possibility is to allow our polyplet fields to decay to SM particles through 
the introduction of higher-dimensional operators which break pseudo-lepton number; in this 
case, we can also relax the requirement that the lightest new particle be neutral. Although 
imposing the pseudo-lepton-number symmetry and then breaking it is somewhat ad hoc, 
such a situation is much more reasonable for more general polyplet scenarios; if n is large, 
the DM stability is ensured by an accidental symmetry, instead of the imposition of pseudo- 
lepton number; in these cases, DM decay would be permitted, but could only occur through 
higher-dimensional operators. Of course, in this case, we evade direct detection constraints, 
but we fail to have a DM candidate. 



V. CONCLUSIONS AND OUTLOOK 

In this work, we have explored a simple "polyplet" extension of the SM which contains 
new fermion SU{2) multiplets and which has the potential to generate a viable DM candi- 
date. This model contains four new Dirac fermions with weak-scale masses and causes no 
difficulties with triangle anomalies. We also find that this simple polyplet model is compat- 
ible with current collider and electroweak constraints and that a DM candidate is possible if 
we are willing to accept the possibility of a heavy {uih > 600 GeV) Higgs boson; addition- 
ally, we expect that solutions where A'^i is a subdominant component of DM can be found 
in the case of a light Higgs boson, although additional fine-tuning may be needed. 

We briefiy comment here on a possible resemblance between the polyplet model and 
supersymmetric models. It may be noted that our n = 2 model contains an SU{2) triplet, 
two doublets, and a singlet; the same is true of the gaugino and Higgsino sectors of the 
Minimal Supersymmetric Standard Model (MSSM). However, several differences exist, the 
most obvious being the lack of sfermion analogues in the polyplet model. Also, all masses 
for the polyplet fermions are Dirac, while those in the MSSM are Majorana. (However, for 
an example of an R-symmetric supersymmetric model-the MRSSM-with Dirac gauginos, 
see lisj.) Lastly, here the Yukawa couphngs are arbitrary parameters; they are not related 
to gauge couplings as in supersymmetric scenarios. 

This work has many possible extensions. Although here we have only considered the 
simplest viable {n = 2) polyplet model, generalizing to higher n would automatically give 
the stabilization of DM, which here was done using a pseudo-lepton-number. Several issues 
would have to be investigated for a specific model with n > 2, such as confirming that 
Yukawa couplings can be chosen such that a neutral member of the polyplet is the lightest 
and calculating the effects of larger multiplets on the oblique parameters. Also, while in 
this work we have insisted that all fields in our polyplet receive their masses from Yukawa 
couplings with the SM Higgs, the phenomenology would likely change substantially if we were 
to allow the introduction of Majorana or dimension-three Dirac mass terms. Additionally, 



One could alternatively introduce a small Majorana mass term for A^i, rendering it an inelastic DM [47 1 
candidate; this would require breaking our pseudo-lepton-number symmetry, but we could still prevent 
A'^i from decaying to SM particles by requiring lepton and baryon number conservation. However, we do 
not pursue this route, as the introduction of an arbitrary mass term not tied to the electroweak scale 
violates the spirit of the polyplet model. 
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it would be interesting to investigate the signatures of polyplets at LHC. We leave these 
explorations for future work. 

Although here we have focused on the relevance of a polyplet model to DM, polyplets are 
interesting in their own right; they are a simple extension of the SM with weak-scale masses 
and reasonably well-defined couplings, which offers the possibility of a rich phenomenology 
accessible at current-day experiments. In polyplet models which do not contain a DM 
candidate, it may be necessary to introduce (possibly higher-dimensional) operators to allow 
the polyplet fields to decay down to SM particles, but these operators may be sufficiently 
suppressed that they have no visible effects at colliders, and the lightest new particle may be 
long-lived on collider time-scales. Additionally, it may be interesting to see if the polyplet 
framework could be used to address other issues, such as neutrino mass. In the n = 2 model 
which we explore here, the SU{2) triplet and both SU(2) doublets can couple to the SM 
leptons at dimension three or four in the absence of pseudo-lepton-number conservation; all 
of these can contribute to neutrino mass. Also, it may be possible to incorporate a polyplet 
containing an SU{2) triplet into a Type-Ill seesaw model. 

In conclusion, fermions placed into larger multiplets of the SM gauge group can be phe- 
nomenologically interesting, yet still undiscovered. The polyplet framework offers a simple 
extension of the SM with implications for dark matter and collider experiments and shows 
promise for future studies. 
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VII. APPENDIX 

Here, we discuss how we can use additional multiplets to loosen the restrictions placed 
on the multiplet quantum numbers by the requirement of anomaly cancellation. We start 
with adding one more copy (or "family") of Xn^.y, fields, i.e.. 



Following the same procedure as those that led to Eqs. ( ITTl) gives the analogous expressions 



(45) 



+ 3(n' 



{n - 
3(n 



(-)(- 



l)A + {n' - 1)A' = 

1 n 1 

- 1)^(2 - yf + 3(- 2 + ^'(^ - 1))(2 - n + - 1)^' 

1 n' 1 

- - Y') + 3(-- + A\n' - !))(- - Y') + {n' - l)^' = (46) 

1 n' 1 

--Y) + f{n- 1)A + (-)(- - Y') + fin' - 1)A' = 0, 



where again A and A' 



can take on the values —1, 0, and 1. 
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We now seek solutions to these equations such that y, y 7^ | to avoid having multiplets 
with zero hypercharge and to avoid having two multiplets with equal but opposite 1^, like 
we saw in the single-family scenario. Let us examine the first of these three equations. The 
solutions to this equation lead us to consider three cases. 

A. Case 1: n = n' 

If n = n', from the first of Eqs. f H6|) . n = n' = 1 or A = —A'. Then, the second and third 
equations yield (| — F) = — (| — Y'). This, however, allows a mass term Xn+i,\-YXn' +i,\-y' 1 
so we discard this option. 

B. Case 2: n 7^ n' , n,n' ^1 

If n 7^ 1, n' 7^ 1 and n 7^ n', from the first equation of fH6l) . A = A' = Q. Then, the 
second and third equation of fH6|) both are fulfilled if(| — F) = — ^(i — Y'). 

C. Case 3: n = 1, n 7^ n' 

For n = 1, n 7^ n', the first of Eqs. ( 146|) gives A' = 0. Then, the second and third 
equations are satisfied if — Y) = —n'{^ — Y'). 

Of course, there is also a case analogous to Case 3 but with the two families interchanged. 

Additionally, we note that, for integer charges, even n implies half-integer Y, while odd 
n implies integer Y; similar relations hold for n' and Y'. For Case 2, this implies that n and 
n' must either both be even or both be odd. For case 3, it implies that n' is odd. 

Of course, more general scenarios are possible; as an example, the SM accomplishes 
anomaly cancellation for a given family using one n = 1 lepton polyplet and three (colored) 
n = 1 quark polyplets. 



[1] A. Knochel and C. Wetterich, Phys. Lett. B 706, 320 (2012) [arXiv:1106.2609' [hep-ph]]. 
[2] R. Foot, H. Lew, R. R. Volkas and G. C. Joshi, Phys. Rev. D 39, 3411 (1989). 
[3] R. Foot, H. Lew and G. C. Joshi, Phys. Rev. D 39, 3402 (1989). 
[4] R. Foot, H. Lew and G. C. Joshi, Phys. Lett. B 212, 67 (1988). 

[5] M. Cirelli, N. Fornengo and A. Strumia, Nucl. Phys. B 753, 178 (2006) [hep-ph/05 12090] . 
[6] R. J. Hill and M. P. Solon, Phys. Lett. B 707, 539 (2012) [arXiv:lllL00T6] [hep-ph]]. 
[7] C. -H. Chen and S. S. C. Law, Phys. Rev. D 85, 055012 (2012) [ arXiv:1111.5 462 [hep-ph]]. 
[8] M. Garny A. Ibarra and S. Vogl, JCAP 1204, 033 (2012) [arXivjlllM155J [hep-ph]]. 
[9] Y. Kajiyama, K. Kannike and M. Raidal, Phys. Rev. D 85, 033008 (2012) [arXiv: 1111.1270) 
[hep-ph]]. 

[10] K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010) and 2011 partial 

update for the 2012 edition. 
[11] R. H. Cyburt, B. D. Fields, K. A. Olive and E. Skillman, Astropart. Phys. 23, 313 (2005) 

,astro-ph/0408033 |. 



20 



[12] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 
[arXiv:1001.4538 [astro-ph.CO]]. 

P. Langacker and G. Steigman, Phys. Rev. D 84, 065040 (2011) [arXiv:1107.3131' [hep-ph]]. 
F. J. Sanchez-Salcedo, E. Martinez- Gomez and J. Magana, JCAP 1002, 031 (2010) 
[ arXiv:1002.3T45] [astro-ph.CO]]. 

C. F. Berger, L. Covi, S. Kraml and F. Palorini, JCAP 0810, 005 (2008) |arXiv:0807.02TT] 
[hep-ph]]. 

P. Batra, B. A. Dobrescu and D. Spivak, J. Math. Phys. 47, 082301 (2006) [h ep-ph/051018T . 
E. Witten, Phys. Lett. B 117, 324 (1982). 

P. Abreu et al. [DELPHI Collaboration], Phys. Lett. B 274, 230 (1992). 

P. J. Fox, R. Harnik, J. Kopp and Y. Tsai, Phys. Rev. D 84, 014028 (2011) [ arXiv:1103.0240] 
[hep-ph]]. 

J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 38, 395 (2005) 
[arXiv:hep-e x/0406019] . 

J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 60, 17 (2009) [arXiv:0901.4486l 

[hep-ex]]. 

Y. Bai, P. J. Fox and R. Harnik, JHEP 1012, 048 (2010) [arXiv: 1005.3797) [hep-ph]]. 
L. M. Carpenter, arXiv:1010.5502 [hep-ph]. 

P. Achard et al. [L3 Collaboration], Phys. Lett. B 580, 37 (2004) [arXiv:hep-ex70310007] . 
T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 85, 011104 (2012) |arXiv:1112. 15771 
[hep-ex]]. 



P. Achard et al. [L3 Collaboration], Phys. Lett. B 517, 75 (2001) [arXiv:hep-ex/0107015' . 
S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710, 26 (2012) [ arXiv: 1202. 1488) 
[hep-ex]]. 

G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 710, 49 (2012) [arXiv:1202.1 408] [hep-ex]]. 

D. Benjamin [CDF and DO and TEVNPH Working Group (The Tevatron New Phenomena 
and Higgs Working Group) Collaborations], arXiv: 11 08. 3331 [hep-ex]. 

J. F. Gunion, ar Xiv:110 5.3965 [hep-ph]. 

L. M. Carpenter, [arXivIlllO, 4895 [hep-ph]. 

M. E. Peskin and T. Takeuchi, Phys. Rev. D 46, 381 (1992). 

Y. Bai, J. Fan and J. L. Hewett, larXiv:1112.T964l [hep-ph]. 

G. D. Kribs, T. Plehn, M. Spannowsky and T. M. P. Tait, Phys. Rev. D 76, 075016 (2007) 
[arXiv:0706.3718 [hep-ph]]. 

M. Baak, M. Goebel, J. Haller, A. Hoecker, D. Ludwig, K. Moenig, M. Schott and J. Stelzer, 

Eur. Phys. J. C 72, 2003 (2012) [arXiv: 1107.0975] [hep-ph]]. 

C. Cheung, M. Papucci and K. M. Zurek. la?Xw:1203.5106l [hep-ph]. 

E. Aprile et al. [XENONIOO Collaboration], Phys. Rev. Lett. 107, 131302 (2011) 
[arXiv: 1104.2549 [astro-ph.CO]]. 

Z. Ahmed et al. [CDMS and EDELWEISS Collaborations], Phys. Rev. D 84, 011102 (2011) 
[arXiv:1105.3377 [astro-ph.CO]]. 

G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) |hep-ph /9506380|. 
P. Agrawal, Z. Chacko, C. Kihc and R. K. Mishra. larXiv:1003. 19121 [hep-ph]. 
G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 180, 747 
(2009) [arXiv:0803.2360 [hep-ph]]. 

K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991). 

M. Ibe, H. Murayama and T. T. Yanagida, Phys. Rev. D 79, 095009 (2009) [arXiv:0812.0072] 



21 



[hep-ph]]. 

[44] M. Ibe, Y. Nakayama, H. Murayama and T. T. Yanagida, JHEP 0904, 087 (2009) 

[arXiv:0902.2914 [hep-ph]]. 
[45] M. S. Dvornikov and A. I. Studenikin, J. Exp. Theor. Phys. 99, 254 (2004) |hep-ph/0411085] . 
[46] J. -F. Fortin and T. M. P. Tait, Phys. Rev. D 85, 063506 (2012) |arXiv: 1103. 32891 [hep-ph]]. 
[47] D. Tucker-Smith and N. Weiner, Phys. Rev. D 64, 043502 (2001) [hep-ph/OlOmS] . 
[48] G. D. Kribs, E. Poppitz and N. Weiner, Phys. Rev. D 78, 055010 (2008) [ arXiv:07 12.2039] 

[hep-ph]]. 



22 



